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Monte Carlo data of the two-dimensional Ising spin glass with bimodal interactions are presented 
with the aim of understanding the low-temperature physics of the model. An analysis of the specific 
heat, spin-glass susceptibility, finite-size correlation length, and the Binder ratio is performed to try 
to verify a recent proposal in which for large system sizes and finite but low temperatures the effective 
critical exponents are identical to the critical exponents of the two-dimensional Ising spin glass 
with Gaussian interactions. Our results show that with present system sizes the recently proposed 
scenario in which the two-dimensional Ising spin glass with bimodally distributed interactions is 
in the same universality class as the model with Gaussian-distributed disorder at low but finite 
temperatures cannot be reliably proven. 
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I. INTRODUCTION 

The two-dimensional Ising spin glassi has been the 
subject of numerous numerical studies^i^i^i^i^ ever since 
its introduction by Edwards and Anderson in 1975fi De- 
spite the fact that this canonical model only orders at 
zero temperature, its popularity can be ascribed mainly 
to its ease of implementation and simplicity. Yet it has 
proven to be extremely difficult to establish if the model 
orders at zero or finite temperature and what the exact 
values of the critical exponents are. Today it is well rec- 
ognized that the model only orders at zero temperature 
In addition, the critical exponents at T = are known 
to good precision^ 

When the interactions between the spins are drawn 
from a Gaussian distribution with zero mean, the system 
has a unique ground state (up to a global symmetry). 
This in turn implies that the critical exponent of the cor- 
relation function has to be zero, i.e., rj — 2 ~ d — 0, 
where d = 2 is the space dimension. In addition, ex- 
tensive zero-temperature domain- wall computations have 
established that the domain-wall stiffness exponent is 
6 = —0.282(2)^ Because at a zero-temperature tran- 
sition the critical exponent of the correlation length v is 
related to 6 via ly — —I/O, it follows that = 3.54(2). 

The ground state of the model with bimodally dis- 
tributed random bonds, on the other hand, is strongly 
degenerated'^" with a finite residual entropy at T = 0. 
Domain-wall stiffness measurements at zero temperature 
show an exponent 6 = (exponential scaling) although 
with significant corrections to scaling up to linear system 
sizes L « 100 under periodic/free boundary conditions. 
Direct measurements of the spin-spin correlation function 
G{r) at T = indicate a critical exponent rj « 0.14.-'^^ An 
extrapolation from finite temperatures^ yielded to high 
precision r/ = 0.138(5). Other estimates of ry have given 
positive values in the range 0.14-Q.40<^iiSii^ii^iii 



Judging from their zero-temperature (critical) proper- 
ties the two versions of the model thus are in two dif- 
ferent universality classes (see Table |n|. However, for 
low but finite temperatures the situation is less clear-cut. 
For the model with bimodal disorder, when a sample is 
in a ground state, turning over one individual spin ei- 
ther leaves the total energy of the sample unchanged or 
increases it by 4J or 8J, where J is the characteristic 
energy of the model (see below). This suggests imme- 
diately that there is an energy gap of AJ between the 
ground state and the first excited state, with A = 4. 
It has been argued however that the true elementary 
excitations are not single spins but more complex ob- 
jects, leading to an effective gap of ~ 2Ji^iii^ii^ Be- 
cause of the gapped excitation spectrum we expect a 
crossover temperature T*{L), which separates a critical 
behavior, in accordance with the aforementioned zero- 
temperature critical point and a critical behavior which 
would resemble a model with continuous interactions. If 
the ground-state degeneracy is given by Nq{L) and the 
excited states have a degeneracy Ni{L), N2{L), . . . , for 
Nq{L) > 7Vi(L)exp(-4J/r) the system will spend al- 
most all of its time in the ground state and equilibrium 
properties will be essentially those of the ground-state 
manifold, whereas for Ni{L) exp{-AJ/T) > Na{L) the 
system will stay in the quantized series of excited states 
and the properties of the system can be expected to 
resemble those of a system with no gap in the energy 
spectrum. Note that this argument is slightly oversim- 
plified as the effects of higher excited states on T*{L) 
are not taken into account, yet we expect their con- 
tributions to be smalliiS, One obtains to lowest order 
T*{L) - 4J/\n[Ni{L)/No{L)]. The ratio Ni{L)/NoiL) 
with a gap 4 J has been estimated in Refs. d and JL^ and 
shows that T*{L) drops gradually as L increases. Our 
data presented below confirm this behavior. 

It has been strongly arguedi^ that in the limit T > 
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T*{L) and with L ^ oo [meaning T*{L) tending to zero 
but never reaching zero] the model with bimodal disorder 
has effective critical exponents identical to the critical ex- 
ponents of the model with Gaussian-distributed disorder, 
so that the two models can be considered as being in the 
same universality class except for the singular behavior of 
the model with bimodally distributed disorder at T = 0. 

We present the results of Monte Carlo simulations of 
the two-dimensional Ising spin glass with bimodally dis- 
tributed disorder on system sizes larger than those used 
m Refs. i and [TtI . Our results show that with current 
algorithms and computer power the data do not pro- 
vide definitive limiting values for the critical exponents 
of the model, although power-law scaling seems plausible 
for finite but nonzero temperatures. Therefore the claim 
that the aforementioned model is in the same universal- 
ity class at finite but nonzero temperatures as the model 
with Gaussian, gap 1/4, or diluted interactionsi^ remains 
to be proven. 

The paper is structured as follows: In Sec.|TT]we intro- 
duce the model, numerical method, and observables, and 
discuss different finite-size scaling relations. In Sec. IIIII 
we summarize previous results on the two-dimensional 
Ising spin glass with bimodally distributed disorder. Re- 
sults on the different critical exponents are presented in 
Sec. IIVI and a finite-size scaling analysis of the data is 
presented in Sec. |Vl 



II. MODEL, OBSERVABLES, AND FINITE-SIZE 
SCALING RELATIONS 



The Hamiltonian of the two-dimensional Ising spin 
glass is given by 



n 



I J Si Sj . 



(1) 



Si ~ ±1 represent Ising spins and the sum is over near- 
est neighbors on a square lattice with periodic boundary 
conditions. The interactions Jij € {^J} (here J = 1) are 
bimodally distributed. For the Monte Carlo simulations 
we use a combination of single-spin flips, exchange Monte 



Carlo updates 



1809 



and rejection-free cluster moves to 



speed up equilibration. Equilibration of the method is 
tested by performing a logarithmic data binning of all 
observables, and we require that the last three bins agree 
within error bars and are independent of the number of 
Monte Carlo sweeps iVs„ccp- The parameters of the sim- 
ulation are listed in Table [H 

The second-moment finite-size correlation 
Iength2a.2i^i23^i25 is given by 
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2 sin(|ki] 



Xsg(O) 

XSG(kmin 
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(2) 



where k,. 



: {2tt/L, 0) is the smallest nonzero wave vec- 
tor, and xsG(k) is the wave- vector-dependent spin-glass 



TABLE I: Parameters of the simulations. A'samp represents 
the number of disorder realizations computed; A'^swcop is the 
total number of Monte Carlo sweeps of the 2Nt replicas for 
a single sample. A'^r is the number of temperatures in the 
exchange Monte Carlo method and Tmin represents the lowest 
temperature simulated. (For L = 128 no data for the specific 
heat has been generated.) 
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susceptibility, 
XSG(k) 



„ik-(Ri-Rj) 



(3) 



In the previous equation [■ • ■ ]av represents a disorder av- 
erage and (• • • ) a thermal average. The finite-size corre- 
lation length is expected to scale as 



(4) 



where i> is the critical exponent for the correlation length. 
This scaling behavior is expected to also be valid for zero- 
temperature transitions when the ground state is not de- 
generate. Since in this work we want to study the ther- 
modynamic limit at finite but nonzero temperatures, we 
postulate that the scaling ansatz in Eq. ^ also holds 
for the model with a bimodal disorder distribution when 
T > 0. 

The standard spin-glass susceptibility XSG = XSG(k — 
0) can also be defined via xsG = ^[('?^)]av, where 



1 ^ 

^ ST^ Qaab 
i=l 



N 



(5) 



In Eq. ([5]) {Sf} and are two copies of the system 

with the same disorder. According to finite-size scaling 
we expect that 



XSG ^ {T ~Tc 



(6) 



and at criticality Xsg{T — T^) ^ L'^^^, where 77 is the 
anomalous dimension exponent of the correlation func- 
tion G{r), 



G(r,r) = [{S.S.+r)' 



(7) 



In addition, we study the dimcnsionlcss Binder ratio^ 
defined via 



g 



mi 



(8) 



In Sec.|V]we plot the Binder ratio as a function of the cor- 
relation length divided by the system size. The method 
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has the advantage that if data for different disorder distri- 
butions he on the same universal curve, the systems are 
in the same universahty class.— Finally, we also compute 
the specific heat of the system)^ 



which is expected to scale as 

Cv - (T-Tc)-". 



(9) 



(10) 



For zero-temperature transitions the critical contribution 
of the specific heat can also be written as Cy{T) ~ T''" 
using zero-transition-temperature scaling relations .-^^^ 
If two systems are in the same universality class they 
share identical values of the critical exponents, as well 
as the values of different observables at criticality [e.g., 
(/(Tc)].— Assuming the power-law behaviors for the dif- 
ferent observables [Eqs. ([4]), ([6]), and (fTO|) ] at low but 
nonzero temperature, we study the values of the effective 
critical exponents. 



III. SUMMARY OF PREVIOUS RESULTS 

The critical properties of the two-dimensional Ising 
spin glass with Gaussian-distributed interactions are 
firmly established from zero-temperature simulations. 
Because the ground state is not degenerate, the corre- 
lation function G{r,T) = 1 for all r at T = 0. 
Therefore, by definition, 77 = 0, ^(T = 0) = 00, and for 
the Binder ratio at zero temperature g{T = 0) = 1. It 
is now well established from domain-wall measurements 
at T = 0, confirmed by size-dependent ground-state en- 
ergy measurements^ that 6 = —0.282(2), hence the ther- 
mal exponent 1/ (whose value is not fixed by the unique 
ground-state condition) is = —1/6 = 3.54(2). Measure- 
ments of v at finite temperatures via Monte Carlo simula- 
tions give consistent estimates;^ yet only if large enough 
system sizes are simulated. The different expected crit- 
ical exponents at zero temperature are summarized in 
Table ni 

For the two-dimensional Ising spin glass with bimodal- 
distributed interactions the situation is, however, much 
more complicated because of the highly degenerate 
ground state, as well as the quantized energy spectrum. 
As noted for instance in Ref. |^, the naive prediction 
for the low-temperature limit specific heat for the two- 
dimensional system of size L and gap A = 4 is 



CviT) 



16 



(rL)2 NoiL) 



(11) 



A priori this low-temperature exponential finite-size be- 
havior should always hold for T <C T*{L). Surpris- 
ingly, there has been a longstanding controversy concern- 
ing the low-temperature behavior of the specific heat in 
the thermodynamic limit. An exponential scaling of the 



TABLE II; Critical exponents for both disorder distributions 
at zero temperature. For the bimodal disorder distribution, 
exponential scaling is expected, i.e., v = 00 (see the main text 
for details). The critical exponent ri in the bimodal case fol- 
lows from Refs. [l^. [^. Isl. [Tol. and [3 S'lid are estimated by 
extrapolating finite-temperature data to T = 0. In the Gaus- 
sian case the estimate for v is from Ref. @. The remaining crit- 
ical exponents can be computed from the zero-temperature 
scaling relations a = —dp and 7 = i^{d — rj), where d = 2 is 
the space dimension. 



Disorder 
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Gaussian 


3.54(2) 





Bimodal 


CO 


0.14 - 0.40 



free energy, and thus correspondingly of all thermody- 
namic quantities, has been first proposed by Wang and 
Swendseni^ They surmised that 



Cy 



-A.J/T 



(12) 



The numerical parameters A = A and P = 2 can be ex- 
pected from the aforementioned arguments regarding the 
gap in the excitation spectrum. In addition, according to 
hyperscaling, the singular part of the free energy scales 
as ^^"^ with d = 2 and so, if Cy scales exponentially, we 
expect that the correlation length ^ scales as 



2J/T 



(13) 



with n = A/2, as predicted first by Saul and Kardari^ 

Wang and Swendsen^ calculated numerically the spe- 
cific heat of the model. Surprisingly, they found A = 2, 
thus suggesting a nontrivial scaling of the free energy. 
However, their measurements were restricted to small 
system sizes and to few disorder realizations, and their 
results indicated strong corrections to scaling. These con- 
clusions stand in contrast to those from work by Saul and 
Kardar,'^ who argue that A = A. In addition, these au- 
thors also estimated n = 2, a behavior which appeared to 
be confirmed independently in work by Houdayer— who 
studied the finite-size scaling of the Binder ratio^^ as well 
as by Katzgraber et al^ who were the first to study the 
finite-size scaling of the finite-size correlation length di- 
rectly via Monte Carlo simulations. Note that an ex- 
ponential scaling of the correlation length (for all A) im- 
plies that the effective critical exponent v is infinite. The 
Wang and Swendsen^ value for A {A = 2) was strongly 
supported by the numerical work of Lukic et al.^ who 
computed the specific heat of the model for intermedi- 
ate system sizes (L < 50), using Pfaffian matrix alge- 
bra techniques. From their analysis they concluded that 
A could be estimated very accurately and that Cy{T) 
tends to the functional form in Eq. (|12p with P = 2 and 
A = 2.02(3),— in agreement with the results of Ref. [j. 
Recently, a new scenario has been proposed in Ref. Il7 : 
While at zero temperature the model with bimodally dis- 
tributed couplings still exhibits exponential scaling with 
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A = 4, at finite but low temperatures in the thermo- 
dynamic Hmit the two-dimensional Ising spin glass with 
bimodal couplings falls into the same universality class 
as the system with Gaussian-distributed disorder. For 
nonzero temperatures and in the thermodynamic limit 
the observables are claimed to display power-law sin- 
gularities with the same critical exponents as the sys- 
tem with Gaussian disorder. In this work we compute 
efi'ective critical exponents with systems larger than in 
Ref. [l3 in the temperature range where power-law scal- 
ing is expected to occurii (0.2 < T < 0.5 for L > 50) 
and show that if the observables can be interpreted to ex- 
hibit power-law scaling, the critical exponents for the sys- 
tem with bimodally distributed couplings cannot be esti- 
mated reliably^'^ from the system sizes studied in Ref. 



IV. RESULTS 

We present numerical data on various observables, 
comparing where appropriate the systems with bimodal 
and Gaussian disorder. We note that for the model with 
bimodal disorder in addition to the crossover temper- 
ature T*{L) due to the energy gap, there is a finite- 
size crossover temperature T^{L) fixed by the condition 
L ^ £,{T). Above T^{L), the observables are close to the 
thermodynamic limit values while below T^{L), XSG a-nd 

tend to be size-limited and thus temperature indepen- 
dent. For all system sizes studied we find T^{L) > T*{L), 
meaning that as T is lowered the size limited condition 
on and xsG sets in well before the effect of the gap. 
The specific heat is approximately size independent down 
to T*, but is strongly affected by the gap: below T*{L) 
the specific heat drops exponentially with decreasing T . 

We now test the hypothesis of Ref. 17 that the effective 
finite-temperature exponents ?7off and exist and are 
identical to the values with Gaussian disorder in d = 2, 
which are ry = and v = 3.54(2), respectively. These 
values imply that the exponent of the specific heat is 
a ~ —7.1 (a — ~dv) and for the susceptibility exponent 
7 ~ 7.1 [7 = v{d - 77)]. 

A. Specific heat 

Figure [T] shows a log-log plot of the specific heat 
Cy{L,T) as a function of T. The numerical results 
are consistent with those of the analogous plot shown 
m Fig. 4 of Ref. Tt', but the present data extend to 
L = 96. The slope in the range 0.2 < T < 0.5 is 
-d\n[Cv{L,T)]/d\n[T] - -4.21(2), i.e., a « -4.21. 
We also perform a point-by-point differential of the data 
for all T (using a second-order midpoint differentia- 
tion combined with a bootstrap analysis to estimate 
the error bars) and thus estimate the effective expo- 
nent aee{L,T) as a function of temperature. Figure [2] 
displays aee{L,T) = -d\n[Cv{L,T)]/d\n[T] as a func- 
tion of T. acf{{L,T) should tend to the thermodynamic 



TT""' — I ' J^' I 




0.1 0.2 0.4 0.6 0.8 1 

T 



FIG. 1: (Color online) Log-log plot of the specific heat Cv 
vs temperature T for several system sizes up to L = 96. For 
0.2 < r < 0.5 the data seem to be approximately independent 
of system size and follow a power-law behavior, i.e., Cy ~ 
T"" with a ^ -4.21(2) (the dashed line with slope 4.21 is 
a guide to the eye). For the system sizes studied, the data 
seem incompatible with a low-temperature, large-size- limiting 
effective exponent a equal to the Gaussian estimate for the 
critical exponent a ~ —7.1 (dotted line with slope 7.1). The 
fluctuations at T < 0.2 can be ascribed to the exponential 
behavior at low enough temperatures. 



critical exponent acB in the limit L — > 00 followed by 
T 0. An extrapolation of the data for T > T*{L), 
where —d\n[Cv{L,T)]/d\n[T] is independent of system 
size L, cannot be performed in a reliable way to test if 
the effective exponent agrees with the expected Gaussian 
value of Qfoff ~ —7.1. 

For 0.2 < T < 0.5 the data seem to saturate (see Fig. [5]) 
although strong fluctuations are present. This "plateau 
region" resembles the behavior predicted by Fiscbi^ who 
argues that there should be a rather broad region in tem- 
perature just above T*{L) for large (but not infinite) L, 
where Cy ~ with x = 5.25(20). This is equivalent to a 
temperature-independent aes = —5.25. Given the uncer- 
tainties in the estimate of the exponent, the agreement 
between the prediction by Fisch and the Monte Carlo 
data presented here is reasonably good. For the system 
sizes studied, which are larger than the ones studied in 
Ref. [TtI . the exponent in this large-L low-T region 
seems to be different from the Gaussian critical expo- 
nent a — ~2v K, —7.1. (In fact, in the Gaussian model 
the low-temperature specific heat is dominated by non- 
critical contributions^ and the true critical behavior is 
not directly visible). Therefore, if both models share the 
same universality class, the system with bimodal disorder 
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FIG. 2: (Color online) Effective exponent acs{L,T) — 
—dln[Cv{L,T)]/d\n[T] as a function of temperature for dif- 
ferent system sizes L. For the system sizes shown, the data 
cannot be extrapolated in a reliable way to OcB = —7.1, the 
value obtained for the model with Gaussian-disordered bonds. 
Therefore it is rather difficult to test if the effective critical 
exponents agree with the known Gaussian values. 

displays huge corrections to scaling and thus simulations 
at considerably larger system sizes would be required to 
prove this beyond any reasonable doubt. 



B. Correlation length and susceptibility 

From the assumed power-law critical behavior of 
the spin-glass susceptibility one can define an effec- 
tive exponent 7eff(i, T) = —d\n[xsG{L, T)]/dln[T] using 
point-by-point differentiation. Figure [3] shows 7efi(i,T) 
against T. The thermodynamic limit domain where the 
effective exponent is size independent can be seen clearly; 
for each L the data peel off the thermodynamic limit 
line at the finite-size limited T^{L). However, even with 
data up to L — 128 there is no reliable way to ex- 
trapolate to the critical value at infinite L and T tend- 
ing to zero. A similar conclusion can be reached for 
jy^ff (L, T) = -dln[C(L, r)]/dln[T], Fig.gl The error bars 
in Figs, m [21 and m have been calculated via a bootstrap 
estimate 1^ 

In Ref. [13 the Caracciolo finite-size scaling 
technique^i^ was used to extrapolate ^{L,T) and 
Xsg{L,T) data for the ±J model towards infinite size, 
in order to extract the exponent 77 from the critical 
scaling relation xsg{L,T) ~ ^(I/,T)^^''. The data are 
interpreted as showing that rj ~ OJ^ 

This technique assumes that there are well-behaved 
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FIG. 3: (Color online) Effective exponent 7cff(I/,T) = 
—dln[xsG{L,T)]/dln[T] as a function of temperature for 
different system sizes L. An extrapolation to the low- 
temperature regime is difficult with system sizes limited to 
L < 128. The horizontal dashed line corresponds to the ex- 
pected Gaussian value 7 « 7.1. An agreement or disagree- 
ment with 7 ~ 7.1 cannot be ruled out. 
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FIG. 4: (Color online) Effective exponent Vb«{L,T) = 
— dln[^(L, T)]/dln[T] as a function of temperature for differ- 
ent system sizes L. An extrapolation to the low-temperature 
regime is difficult with system sizes limited to L < 128. Any 
extrapolation to the low-T behavior where v « 3.54 would be 
difficult to perform with the current data. 
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FIG. 5: (Color online) i{2L,T)/^{L,T) against ^{L,T)/L. 
If corrections to scaling were small, all data should lie on an 
universal curve. This is a prerequisite to perform an extrap- 
olation to the bulk regime. Clearly, the data do not scale 
well. 



scaling functions for ^{L,T) and XsGiL,T), such that, 
in particular, ^(2L, r)/^(L, T) is a unique function of 
£,{L,T)/L. Otherwise the Caracciolo scaling procedure 
cannot be applied as systematic errors would be intro- 
duced. Figure M shows data for ^{2L,T)/^{L,T) vs 
£,{L,T)/L with L between 16 and 48. It can be seen 
that the finite-size corrections to scaling are strong be- 
cause the curves for different L do not superimpose, and 
also because the ratios £,{2L, T)/^{L, T) increase beyond 
a value of 2 at low temperatures (in the absence of cor- 
rections 2 is the strict Tc limit of the ratio) . This means 
that the Caracciolo procedure (or any similar protocol 
such as the one by Kim-i used by Katzgraber et al^ for 
the model with Gaussian disorder) has to be performed 
with considerable care for bimodally distributed disorder 
in two space dimensions. In Fig. [HI we show a log-log plot 
of the correlation length as a function of temperature to- 
gether with an extrapolation using the methods of Kim^^ 
as well as Palassini and Caracioloi22i2^ Both extrapola- 
tions agree very well. While the extrapolated data seem 
to follow a power-law behavior with v w 3.45, the extrap- 
olation method is not reliable due to strong corrections 
to scaling (see Fig. [S])- Until data on much larger system 
sizes become available it does not seem plausible to give a 
reliable account of the large-L, low-T limiting functional 
behavior of xsciL, T) or ^(L, T). 

Alternatively, with no extrapolation, one can define an 
effective exponent 



FIG. 6: (Color online) Log-log plot of the finite-size corre- 
lation length ^{L,T) as a function of temperature for sev- 
eral system sizes. The solid orange circles represent extrapo- 
lated data to the thermodynamic limit from L = 96 using the 
method of Kirti^i and the solid pink triangles represent the 
data extrapolated to L = oo using the method of Palassini 
and Caracciolo (PC).— The (extrapolated) data seem to fol- 
low a power-law behavior with f ^ 3.45, which is close to the 
value of the critical exponent for Gaussian-distributed disor- 
der, f = 3.54(2). The dashed line is a guide to the eye. 



2-r/eff(i,T) 



d\n[xsG{L,T)] 
dlnK(L,r)] 



(14) 



via differentiation of the data. For the available system 
sizes rics{L,T) is always greater than ^ 0.2 (see Fig. [T]). 
A reliable extrapolation to infinite L and T tending to 
zero would again require data of much larger system sizes. 

There have been numerous estimates of rj at zero tem- 
perature for the two-dimensional Ising spin glass with 
bimodally distributed disorder, including direct measure- 
ments of the correlation function G{r) by combinatorial 
or Monte Carlo methodS )^°i^^ and there is a general con- 
sensus that 77 > 0.15. Indeed it can be noted that the 
finite-T raw G{r) data in McMillan's^^ Fig. 1 are by in- 
spection incompatible with = 0. 

From the definition of 77 through G{r), when there is 
a degenerate ground state at T = 0, the time average 
spin-spin correlation function must decay with increasing 
distance r implying a positive-definite value for rj. A 
fortiori at any T slightly above zero one would expect 
G{r,T) < G{r,0) for all r except in quite exceptional 
cases. Otherwise, a limiting rj^g = [meaning G(r, T) = 
1 for all r at T close to zero] appears to be ruled out from 
basic physical principles. 
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FIG. 7: (Color online) Effective exponent 2 — ricti{L,T) — 
dln[xsG{L,T)]/dln[^{L,T)] as a function of temperature for 
different system sizes L. For all system sizes and temper- 
atures studied rj^s is always greater than 0.2, although an 
extrapolation to rjctt — cannot be ruled out. 



C. Gaussian disorder 



We have also computed the effective critical exponents 
for the correlation length [foff(£,T)], the specific heat 
[aes{L, T)], the susceptibihty [7eff(L, T)], and correlation 
function [r^off (i, T)] for the two-dimensional Ising spin 
glass with Gaussian-distributed disorder in order to test 
corrections to scaling in that model. In this case 



P(J.,)~e-^i/2'^ (15) 



in Eq. ([T]). The results are qualitatively similar to the 
results found for the model with bimodally distributed 
disorder, but the data extend to lower temperatures thus 
making an extrapolation to zero temperature slightly 
more reliable. Still, without the knowledge of the zero- 
temperature estimate of the stiffness exponent 9, v ~ 
— 1/9 could not be determined to such high precision. In 
Fig. [5] we illustrate this behavior with data for T^cff and 
I'eft as a function of temperature. In Fig. [5] we compare 
the effective critical exponent for Gaussian and bi- 
modal disorder for L — 128. While the extrapolation in 
the Gaussian case can be done easily to T = 0, since the 
zero-temperature limit is well known, this is difficult for 
the bimodal case. 




FIG. 8: (Color online) Effective exponent 2 — ricti{L,T) = 
dln[xsG{L,T)]/d\n[(,{L,T)] as a function of temperature for 
different system sizes L for Gaussian disorder. The data ex- 
trapolate well to 77 « 0. Note that for L = 128 the data 
for T > 0.9 have been dropped due to strong fluctuations. 
Inset: Effective exponent VeS as a function of temperature 
for different system sizes L. The dashed line corresponds to 
the zero-temperature estimate from the stiffness exponent, 
u = 3.54(2). 



V. UNIVERSALITY AND FINITE-SIZE 
SCALING 

The critical values of the correlation length divided by 
the system size ^l(Tc)/L [Eq. ^] and the Binder ratio 
g{Tc) [Eq. ([H])] are characteristic of the universality class 
of a continuous transition. These are linked to 77 at crit- 
icality because they represent various ratios of integrals 
with G(r) ^ Pqj. example, in a strip geometry 

for two space dimensions, S,l/L = l/iirr]) at criticality^^^ 
The correlation length divided by system size of the two- 
dimensional spin glass with Gaussian-distributed disor- 
der diverges for T ^ and the Binder ratio tends to 
1. If the two-dimensional Ising spin glass with bimodal 
interactions lies in the same universality class as the sys- 
tem with Gaussian disorder, then identical values for 
these parameters at criticality should be observed. The 
Binder ratio values g{L,T) become temperature indepen- 
dent within the error bars for T < T*{L) providing esti- 
mates of the zero-temperature values g{L,0). There are 
corrections to scaling but the series of points appear to 
tend to a large-L limit which is significantly less than 
unity. 

For any continuous transition, at large L and T ap- 
proaching Tc the Binder ratio is a nontrivial function of 
the variable ^l/^-^'^'^^'— A plot of g{L, T) vs Cl(L, T)/L 
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FIG. 9: (Color online) Efltective critical exponent l/7off for 
both Gaussian (open symbols) and bimodal disorder (full 
symbols) as a function of temperature T. While in the Gaus- 
sian case the T = limit is known (7"^ ~ 0.14, black square) 
and thus the extrapolation from the finite-temperature data 
can be performed via a simple third-order polynomial (dashed 
line), the data for the bimodal case can be extrapolated to 
any arbitrary value including l/7cff = 0, which corresponds 
to exponential scaling. 



FIG. 10: (Color online) Binder ratio g{L, T) as a function of 
the finite-size correlation length divided by the system size, 
S,l{L,T) / L. While the data for Gaussian disorder seem to 
extrapolate to g{L,0) — 1 and ^l{L,T)/ L — 00, the data 
for bimodal disorder seem to stop at g{L,0) ~ 0.92 and 
£,l{L,T)/ L « 0.91 (marked by a boxed cross). 



VI. CONCLUSION 



(Fig.[TU| shows a remarkable behavior: aU points for both 
bimodal and Gaussian disorder are on a unique curve. 
This is particularly striking as the data span both the 
regions T > T^{L) and T < T^{L). There is, however, 
a qualitative difference between the Gaussian data and 
the bimodal data. For the former at each L the data 
points extend to the same zero-temperature end point 
\^{L,T = 0)/L = 00 and g{L,T = 0) = 1, while for the 
system with bimodal disorder the end points for differ- 
ent L seem to cluster and not grow beyond ^l{L,T — 
0)/L ~ 0.91(2) and g{L,T = 0) ~ 0.92(2). In fact, 
the data do grow slightly, but the growth rate is within 
statistical error bars. This point can be interpreted as 
the zero-temperature critical parameters for this model. 
Thus, while the scaling functions agree, at zero tempera- 
ture both models seem not to be in the same universality 
class. An alternative explanation could be that the model 
with bimodal disorder is "marginal," i.e., the endpoint 
might approach g{L,T = 0) = 1 logarithmically slow. 
In the bulk regime, which corresponds to the lower left 
corner of Fig. [TOl data for g{L,T) and £,l{L, T)/L agree 
and thus suggest that both models might share a common 
finite-temperature universality class. Note that this can- 
not be inferred from studying the critical exponents due 
to large corrections to scaling, as shown in Sec. IIVI 



We have carried out a careful numerical study of the 
two-dimensional Ising spin-glass model with bimodally- 
distributed disorder at finite temperatures with system 
sizes up to L = 128, measuring the specific heat, finite- 
size correlation length, the spin-glass susceptibility, and 
the Binder ratio. The data neither confirm nor disprove 
the hypothesisi^ that in the large-L, low-but-nonzero-T 
limit the effective critical exponents of the system are 
identical to the known critical exponents of the two- 
dimensional Ising spin glass with Gaussian-distributed 
disorder, although the data seem to follow a power-law 
behavior at finite T (see, for example. Fig. [J). Since the 
evidence for the critical exponents of the two-dimensional 
Ising spin glass with bimodally distributed disorder be- 
ing identical to those of the pure Gaussian case is weak, 
a further critical analysis of the estimates for the other 
distributions studied in Ref. [H might be opportune. A 
plot of the Binder ratio against the finite-size correlation 
length (Fig. [TO]), which shows very small corrections to 
scaling, only suggests that in the bulk regime and for fi- 
nite temperatures both models might share a common 
universality class. Corrections to scaling are extremely 
large in the bimodal case. Therefore our main result is 
that simulations with larger system sizes and lower tem- 
peratures are needed to conclusively compute the criti- 
cal exponents to prove the scenario proposed in Ref. [l3 
in which both models with Gaussian and bimodally dis- 
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tributed disorder are in different universality classes at 
T = 0, where tfie degeneracy of tfie ground state plays 
a key role, yet the models share the same universality 
class at finite nonzero temperatures and very large sys- 
tem sizes. 



ful discussions, and T. Jorg for pointing out the useful- 
ness of plotting the data in the way shown in Fig.[TOl The 
simulations were performed on the Hreidar and Gonzales 
clusters at ETH Zurich. 
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